Reconstructing the interaction between the dark matter and holographic dark energy 
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We reconstruct the interaction rate between the dark matter and the holographic dark energy 
with the parameterized equation of states and the future event horizon as the infrared cut-off length. 
It is shown that the observational constraints from the 192 SNIa and BAO measurement permit the 
negative interaction in the wide region. Moreover, the usual phenomenological descriptions can not 
describe the reconstructed interaction well for many cases. The other possible interaction is also 
discussed. 

PACS numbers: 95.36.+X, 98.80.-k, 98.80.Cq 



I. INTRODUCTION 



In the modern cosmology, a 'dark energy' (DE) with negative pressure is suggested to be responsible for the current 
acceleration of the universe. The simplest candidate of DE is the cosmological constant, which does nicely well at the 
pragmatic observational level, but entails the serious theoretical difficultly: the cosmological constant problem and 
the coincidence problem. Explanations of DE have been sought within a wide range of physical phenomena, including 
some exotic fields, modified gravity theories, and so on -see [l| and references therein. Among the most recent generic 
proposals, the model inspired by the holographic idea 0, that the quantum zero-point energy of a system cannot 
exceed the mass of a black hole with same size, has been put forward to explain the DE [1, 0, IE]- This DE density 
can be determined in terms of the horizon radius of the universe, corresponding to relate the UV cutoff of a system to 
its IR cutoff in the quantum field theory. There are usually three choices for the horizon radius supposed to provide 
the IR cutoff, with different degrees of success, namely the Hubble horizon, the particle horizon, and the future event 
horizon. The event horizon may be better, since in this case the DE can drive the present accelerated expansion and 
the coincidence problem can be resolved by assuming an appropriate number of e- folding of inflation 0, [a, 0] . 

Most discussions on DE models rely on the fact that both dark matter (DM) and DE only couple gravitationally. 
However, given their unknown nature and the symmetry that would impose a vanishing interaction is still to be 
discovered, an entirely independent behavior between dark sectors is very special. Moreover, since DE must be 
accreted by massive compact objects like black holes and neutron stars, in a cosmological context the energy transfer 
from DE to DM may be small but must be non-vanishing. The interaction hypothesis was first introduced by 
Wetterich Q to discuss the cosmological constant problem in the light of dilatation symmetry and its anomaly. Then 
cosmological consequences of a scalar field coupled to the matter were studied in 0] . It was found that the coupling 
quintessence models may give the scaling attractors providing an accelerated expansion at the present time and 
alleviate the coincidence problem [l(| ■ The interaction also appears in the context of modified gravity models [ll[ . 
More possibility that DE and DM can interact has been studied in [12, 0, [3 HE EE E3 • Confronted to cosmological 
data, it was found that an appropriate interaction can influence the perturbation dynamics, and the lowest multipoles 
of the CMB spectrum [IE and could be inferred from the expansion history of the Universe, as manifested in 
the supernova data together with CMB and large-scale structure [20, HJ, [22| • In addition, it was suggested that the 
dynamical equilibrium of collapsed structures would be affected by the coupling of DE to DM [23|, [24[ . The interaction 
was first connected to holography by Horvat [25| who argued that scaling of the cosmological constant stemming from 
the zero-point energy in quantum field theory possibly implies a non- vanishing coupling of the cosmological constant 
with DM. In the holographic DE model with the Hubble horizon as the IR cutoff the interaction can be available to 
derive the present accelerated expansion and alleviate the coincidence problem [261 ] . In the interacting model with the 
event horizon as the IR cutoff, it was shown that the equation of state (EoS) of DE can accommodate the dynamically 
evolving behavior of crossing the phantom divide [27j , which suggested by recent most observational probes [28[ . 

Although the interaction is important in studying the physics of DE, it will not be possible to derive the precise 
form of the interaction from first principles unless the nature of both dark sectors was known. Usually, the coupling 
is determined from phenomenological requirements [lO, 26] . In view of the continuous equations of DE density pa and 
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DM density p m , the coupling Q must be a function of densities multiplied by a quantity with units of inverse of time, 
which has an obvious choice as Hubble time H . Thus, one may write the coupling as 

Q = Q(Hp m ,H Pd ), (1) 

which leads Q ~ X m Hp m + XdHpd from the first order terms in the power law expansion. Assuming that the ratio 
r = Pm/pd might be piecewise constant, the linear parameters are usually set to A m ~ A^ and even A m ~ or A^ — 
for simplicity. Considering the couplings are terms in the Lagrangian which mix both DE and DM, one may further 
suppose that they could be parameterized by some product of the densities of DE and DM, such as the simplest 
Q ~ \pmPd 0- Besides these phenomenological descriptions, various proposals at the fundamental level have been 
tried to account for the coupling, including the dependence of the matter field on the scalar field [29[ or expressing 
the cosmological constant as a function of the trace of the energy-momentum tensor [301 ] - Recently, an interesting 
thermo dynamical description of interaction between holographic DE and DM has been proposed in [3l|, where it 
was assumed that in the absence of the coupling the DE and DM remain in separate thermal equilibrium, then a 
small interaction can be viewed as a stable thermal fluctuation that brings a logarithmic correction to the equilibrium 
entropy of DE and DM. Other specific coupling which was assumed from the outset can be found in pjl, [32. [33j. 

The main aim of this work is to reconstruct the coupling using the recent DE probes (the Baryon Acoustic Oscillation 
(BAO) measurement at z = 0.35 from the Sloan Digital Sky Survey [34| and the re-compiled 192 Type la Supernovae 
(Snla) samples [33, consisting 60 points from ESSENCE ("Equation of State: Supernovae trace Cosmic Expansion") 
supernova survey [3(|. 57 points from Supernovae Legacy Survey [37|, 45 points nearby Supernovae (38|, and 30 points 
detected by the Hubble Space Telescope [39|). We will focus on the holographic DE model and choose the future 
event horizon as the IR cutoff. 

The model with the Hubble horizon has been studied recently in [4(| and it was found that the reconstructed 
interaction is always positive in la region. This seems to corroborate the recent argument that the negative interaction 
violates the second law of thermodynamics which inquires the energy transfer from DE to DM rather than otherwise 
[4lj . However, it should be noticed that there are some problems in the thermodynamics of DE, such as the negative 
entropy [42j], and the generalized thermodynamical second law indeed breaks down in the universe with phantom- 
dominated DE [43|, [44| . These results suggest that one should consider the thermodynamical properties of DE with 
wide possibilities. Hence, it is interesting to see whether the positive interaction is a robust result for other models, 
such as the present model with the IR cutoff as the future event horizon. 

Considering the time varying DE gives a better fit than a cosmological constant and in particular most of the 
observational probes indeed mildly favor d yna mical DE crossing the phantom divide at z <~ 0.2 (28|, we will employ 
two commonly used parameterizations pll. |45|. l46l. |47|| . namely 

2: 

w A = w + wx(l - a) = w + wi — , (2) 

r 1 z 

which has been used in [40j, and 

wb = Wo + wi(l — a)a = wq + w\ ^. (3) 

{l + zf 

It should be noticed that the different parameterizations are beneficial to control some amount of parameterization 
dependence. After reconstructing the interaction, we will further compare it with the usual phenomenological models 
and the recent thermodynamical description. 



II. RECONSTRUCTION 



Let us begin with the Fricdmann equations 



H 2 = ~ ( Pm + p d ) , (4) 



H = -~(p m + Pd+Pm+Pd), (5) 

where we have normalized 8nG = 1 for conventions. The total energy density p — p m + pd satisfies a conservation law. 
However, since we consider the interaction between DE and DM, p m and pd do not satisfy independent conservation 
laws, they instead satisfy two continuous equations 



Pm + 3Hp m = Q, 



(6) 
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p d + 3H(l + w d )p d = -Q, (7) 

where w d is the EoS of DE, and Q denotes the interaction term. Without loss of generality, we will write the interaction 
as Q = p d T, where T is an unknown function. 

Using the the ratio of energy densities r = p m / Pd, we have 

H 3 1 + w d + r ^ 



from Eq. ([5]). From Eq. ([6]), we have 



H 2 2 l + r 



r = —r— — 3Hr + F, 



which can be recast as 

f = (1 + r)r + 3Hrw d (9) 
from Eq. l[7]l. Eliminating the T in above two equations, we obtain 

-rr 3Hrw d f pa + ^ 



l + r p d 

Until now, we have not specified the concrete DE density. We will focus on the holographic DE model. Followed 
Q by choosing the future event horizon 

/>oo J 

R *=°L it* ("> 

as the IR cutoff, the holographic DE density is p d — 3c 2 i?^ 2 , where c 2 is a constant and the Planck mass has been 
taken as unit. The most possible theoretical value of c is one [El, IB], indicating that the total energy from DE must be 
determined by the Schwarzschild relation. Taking the derivative with respect to t, the evolution of the horizon can 
be determined by 

Re = HR E - 1. 



Defining fl d = p d /3H 2 , we have il d = 1/(1 + r) and Re = cy/1 + r/H. Thus, Eq. (fT0|) can be recast as 

-rr - 3Hrw d _ -2H (c-y/l + r - l) 



l + r cy/1 + r 



3Hr. (12) 



Obviously, r is not a constant in general. This is different with the case in [40] where r is a constant since the IR 
cutoff was chosen as being the Hubble scale. Replacing the time t as the redshift z = (1 — a) /a, we can rewrite Eq. 
p]) as 

-(1 + z)rr' + 3rw d _ 2 (cyT +T - l) _ ^ 



l + r cy/1 + r 

where the prime denotes the derivative with respect to z. Similarly, Eq. {8j reads 

H' 3 l + w d + r 

IT ~ 2(1 + 0) l + r 



(14) 



It is interesting to find that Eqs. (| 13|) and (fi"4"|) determine the evolvement of r and iJ, if we know the EoS w d . In the 
normal interacting DE model, one assumes the explicit interaction form to give out w d . For our aim, we will use the 
two commonly used parameterizations of w d to determine the dynamics of our model, and reconstruct the interaction 
rate T/3H 

T (1 + z)r' + 3rw d 

3H ~ 3 (l + r) 
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FIG. 1: The deceleration parameter q with the best-fit values of equation of states w^. The red, green, and blue lines denote 
the models without interaction, with interaction and Wd = wa, and with interaction and Wd = wb, respectively. For each 
model, there are five lines from top to down at early taking different c from 0.6 to 1.4. 

from the recent observational datasets. 

We will use the re-compiled 192 Snla samples (0 < z < 1.8) combined with the recent BAO measurement from 
SDSS to reconstruct the interaction rate. There are other DE observational probes, including the three-year WMAP 
CMB shift parameter, the X-ray gas mass fraction in clusters, the linear growth rate of perturbations at z = 0.15 as 
obtained from the 2dF galaxy redshift survey, and the look back age data. However, since the parameterizations of 
Wd Eqs. @ and ([3]) are motivated to accommodate the dynamically evolving behavior of crossing the phantom divide 
at recent epoch, and our model has not included the radiation and the baryonic matter which may be important in 
the early, we will not use the WMAP CMB shift parameter which focuses on the high redshift region. Besides, for 
simplicity, we do not adopt other probes of DE which have large relative errors compared with Snla, CMB, and BAO 
probes [48|. 

As usually, we will fix DM density parameter as f^om = 0.3 or 0.25 to include the best-fit value of Qom = 0-27 
from 5-year WMAP data. In general, they are sufficiently representative. Moreover, it is convenient to compare our 
reconstructed interaction to the interaction reconstructed in (40| where these two DM density parameters are used. 
We will consider five indicative different values of the constant c near one. The best-fit values with la error bars 
for the parameters wq and w\ are given in TABLE [JJ [Til One can find that the phantom divide crossing in recent 





c 


0.6 


0.8 


1.0 


1.2 


1.4 


WA 


wo 


-1.19±0.18 


-1.16T0.16 


-1.14±0.16 


-1.13i0.15 


-1.12i0.14 


Wl 


0.60L1.26 


1.09L1.07 


1.35i0.97 


1.50i0.90 


1.60i0.85 


wb 


wo 


-1.21±0.24 


-1.21T0.22 


-1.21i0.20 


-1.21i0.19 


-1.20i0.19 


Wl 


0.99L2.11 


1.83il.83 


2.29il.66 


2.58il.55 


2.78il.47 



TABLE I: The best fit values with la error bars for the parameters of Wd with DM density parameter Qom = 0.25. 





C 


0.6 


0.8 


1.0 


1.2 


1.4 


WA 


w 


-1.03±0.20 


-1.03L0.18 


-1.03L0.16 


-1.03±0.16 


-1.02±0.15 


Wl 


0.48L1.28 


0.48L1.09 


0.48L0.98 


0.70L0.92 


0.85L0.87 


WB 


wo 


-0.99±0.26 


-1.02L0.23 


-1.04L0.22 


-1.05±0.21 


-1.05±0.20 


Wl 


-1.14±2.20 


0.03L1.90 


0.68L1.73 


1.09L1.62 


1.37L1.54 



TABLE II: The best fit values with la error bars for the parameters of Wd with DM density parameter Qom = 0.3. 

epoch is always permitted in la region. The deceleration parameter with these best-fit equation of states is plotted 
in FIG. [TJ where the acceleration (super-acceleration, at most cases) in recent epoch is achieved. It is interesting to 
see that current deceleration parameter in interacting model is almost not affected by the DE parameter c, according 
to the almost same parameters wq in TABLE [TJ [TTJ We also show the evolution of p m compared its non-interacting 
version p^ m in FIG. O The ratios change strongly in recent epoch which implies the interaction plays important role 
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FIG. 2: The ratios Pm/Pm f° r different fiom, Wj, and c. In each panel, c increases from top to bottom. 
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FIG. 3: The ratios Pd/pd f° r different fiom, %, and c. In each panel, c increases from top to bottom. 

in that time. The correspnding plots for the dark energy component are given in FIG. [3l Since the DE density is 
always positive, this suggests that the interaction can be negative but can not be too negative when Wd > —1, which 
is possible at least in the past. 

With these best-fit values, we reconstruct the interaction rates, see FIG. 0] and FIG. [5] This is one of the main 
result of this paper. It should be stressed that the negative interaction is permitted in the wide region of FIG. Q] and 
FIG. [SJ This is contrast to the result in (40| where the reconstructed interaction is alwa ys p ositive, and the argument 
that the second law of thermodynamics which imposes energy transfer from DE to DM [4l| . If the second law is true, 
the permitted region of the interaction rates must be reduced, as done in 49]. However, for generality and considering 
the thermodynamical properties of DE is not clear, we will not restrict the interaction to be positive. 

Moreover, the effective EoS u> e // = Wd + jjj can be obtained since we have known Wd and reconstructed the T/3H. 
From TABLE miTTl and FIG.IUEl one can know that the effective EoS may be bigger or smaller than —1 for different 
parameters flo m and c. If one requires an effective phantom- like DE, the presence of a coupling makes the parameters 
have more possible values than the case of absence of a coupling, where the EoS is determined by Wd = — | — ■^\ / Tld 
which imposes c < y/Tld for Wd < —1. Moreover, it is possible to have effective EoS w e ff smaller than —1 but Wd 
bigger than —1 when the interaction is negative. This suggests that the negative interaction is interesting because 
DE with Wd > —1 is easily accepted (The scalar field with Wd < —1 will break the zero energy condition.) and the 
effective EoS of DE w e ff < — 1 can fit the cosmological data better. 

III. COMPARISON 

In the following, we will compare the reconstructed interaction rate with other descriptions. Let us begin with 
the usual phenomenological descriptions. Usually there are four different choices of Qi (i — 1,2,3,4), which can be 
expressed as 3XHpd, 3\Hp m , 3\H(p m +pd), and Xp m Pd, respectively. To compare them with the reconstructed 
interaction rate T/3H, we define four interaction rates Ri = Qi/ (3Hpd), which can be determined by r and H (which 
have been solved from Eqs. fTB")) and fLi]) ). namely Ri = A, R 2 = Ar, R 3 = A(l + r), and R 4 — XHr/ (1 + r). Observing 
FIG. [H the indicative curves of the usual interaction rates Ri (i = 1,2,3,4) show that in many cases they are not 
favored in la region. For example, see the case of c = 1.2, where the interaction rates R2 (green line) and R3 (blue 
line) are permitted in la region, but the other interaction rates R\ (red line) and R4 (yellow line) are not permitted 
whether A increases or decreases. For the case c = 0.6, the reconstructed interaction rate changes at recent epoch 
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FIG. 4: Interaction rates for flom = 0.3 with respect to redshift z. Above and below panels denote Wd — wa and wb, 
respectively. The black curves and the grey region between them indicate the best fit curve of the reconstructed interaction 
rate T/3H and the region in la confidence level. The black dashing lines indicate T/3H = 0. The red, green, blue, and yellow 
curves indicate the Ri from i = 1 to i = 4 with four constants A. For example, we give A = 0.13,0.11,0.07,0.2 in turn for 
Wd = wa and c = 1. The blue dashing lines indicative the interaction rate R$. The red dashing lines in the case with c = 0.6 
indicative the interaction rate R$ with two nonzero constants. For the above panel, they are Ai = 0.17, A2 = —0.16; For the 
below panel, they are Ai = 0.11, A2 = —0.15. 




0.2 0.5 0.8 1.1 1.4 1.7 0.2 0.5 0.8 l.I 1.4 1.7 0.2 0.5 0.8 1.1 14 1.7 0.2 0.5 0.8 1.1 1.4 1.7 0.2 0.5 0.8 1.1 1.4 1.7 




0.2 0.5 0.8 1.1 1.4 1.7 0.2 0.5 0.8 1.1 1.4 1.7 0.2 0.5 0.8 1.1 1.4 1.7 0.2 0.5 0.8 1.1 1.4 1.7 0.2 0.5 0.8 1.1 1.4 1.7 



FIG. 5: Interaction rates for Qom = 0.25. The blue dashing lines indicative the interaction rate R$. The red dashing lines 
indicative the interaction rate R§ with some constants Ai (« = !,-•■ ,5). 
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from positive to negative. This tendency is more general in FIG. [5] with small DM parameter flom — 0.25. Thus, one 
may suspect that the four phcnomcnological descriptions of the interaction which have the definitive sign may be not 
suitable. 



We will also consider the recent thermodynamical description of the interaction. Defining Q§ — 3Hb 2 (p n 
the thermodynamical interaction rate b 2 can be determined by following equations (see [3l| in detail): 



2il 



3/2 



3c 



m/c-i 



H 2 



Kin 



yi-kc 2 c /^rr d (c/^r d - 1) 



-fiS)'(i+*), 



Pd) 



(«d - 1) 



{Q d - 1) = -3b 2 



~(l + z)H' 



- 1 



H c 2n d ' 

where the zero superscript of tt d and H° indicates absence of interaction. We can relate b 2 to a new interaction rate 

Q5 _ b 2 (p rn + p d ) 



R* = 



SH Pd 



Pd 



0L 

n d 



In FIG. 2] and FIG. [5] (blue dashing lines), one can find that R5 is almost not favored in la region (Note we have 
not shown the cases with c < 1 where b 2 has not the real solution and the parameterization of w d is not needed to 
determine i?s since the EoS in this model is determined by the b 2 .). 

It is interesting to ask whether the reconstructed interaction can be described well by a more general interaction 
form. A natural candidate is to expand the phcnomcnological description (TTJ up to the second order, namely 



Qe = 3XiHp m + 3X 2 Hp d + X 3 H 2 p mPd + MB 2 p 2 m + X 5 H 2 p 2 d , 



with the interaction rate 



Re 



Qe 



3H Pc 



H 3 r H 3 r 2 H 3 

= Air + A 2 + A 3 - h A 4 - h A 5 - 



1 



1 



T +r 



We find that it indeed works well, see the red dashing line in FIG. H] and FIG. [5] However, this parameterization 
only has theoretical interest since it contains too many parameters. Moreover, numerical calculations prove that this 
parameterization can not be reduced to include only the two first order terms, which has been studied recently in 
[49|, if we need it being permitted in la region for all cases. 



IV. SUMMARY 



We have reconstructed the interaction term between the holographic DE and DM, using the re-compiled 192 Snla 
samples combined with the recent BAO measurement. The DE parameter c is assumed near one. The two common 
used parameterizations of w d are considered in la region. It is found that the present accelerated expansion of universe 
is achieved and the phantom behavior of DE is permitted. We illustrate that the negative interaction is permitted 
in the wide region. Hence, contrast to the DE model studied in [4(J, where the reconstructed interaction is always 
positive in la region, we can not obtain the favor from the DM observation for the recent argued thermodynamical 
second law which imposes the energy transfer from DE to DM. This suggest us to keep wide possibilities of the 
thermodynamical properties of DE. 

We show that the four usual phenomenological descriptions can not describe the reconstructed interaction well for 
many cases. Specially for the case with small DM parameter 57om = 0.25, the interaction rate has a trend to change 
sign at recent epoch. This is at variance with the usual phenomenological interacting terms which have the definitive 
sign. We further illustrate that the recent thermodynamical description of the interaction is not favored. Our work 
stimulates one to seek a more suitable interaction. 

It would be interesting to confront our model to more observations, such as CMB angular power and large scale 
structure, which will further constrain the interaction rate and make clear if the negative interaction and phantom DE 
are still to be permitted. Another interesting work is to consider the perturbation evolution in our model. Recently, 
it has been found that some types of interaction will lead to instability under the curvature perturbation [Hoi . l5l| . It 
may provide some restrictions on the form of interaction. We will work on these directions in the future. 
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